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ABSTRACT 

We apply a new algorithm, called the Unbiased Minimal Variance (here- 
after UMV) estimator, to reconstruct the cosmic density and peculiar velocity 
fields in our local universe from the SEcat catalog of peculiar velocities com- 
prising both early (ENEAR) and late type (SFI) galaxies. The reconstructed 
fields are compared with those predicted from the IRAS PSCz galaxy redshift 
survey to constrain the value of /3 = il^f /b, where Q m and b are the mass den- 
sity and the bias parameters. The comparison of the density and velocity fields 
is carried out within the same methodological framework, and leads, for the 
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and 13 = 0.51±0.06, 



first time, to consistent values of (3, yielding /3 
respectively. 

We find that the distribution of the density and velocity residuals, relative 
to their respective errors, is consistent with a Gaussian distribution with a ~ 1, 
indicating that the density field predicted from the PSCz is an acceptable fit 
to that deduced from the peculiar velocities of the SEcat galaxies. 
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1 INTRODUCTION 

In the gravitational instability scenario {e.g., Peebles 
1980), mass density fluctuations and peculiar veloci- 
ties evolve in an expanding universe under the effect 
of gravity. If density fluctuations are small, linear the- 
ory is valid and a simple relation exists between peculiar 
velocities, v, and mass density contrast, S m : 

V • v = -n° m 6 6 m , (1) 

where Q m is the mass density parameter. Equation (jjj) 
shows why peculiar motions are so important in cos- 
mology: they provide a direct probe of the mass density 
distribution in the universe. The mass density fluctua- 
tion field, 8 m , can be deduced from the galaxy observed 
density contrast, <5 9 , assuming a relation (bias) between 
the distribution of galaxies and that of the underlying 
density. The simplest relation suggested in the litera- 
ture is that of linear bias, namely S g — b8 m , where b 
is the linear bias parameter for a given population of 
mass tracers. This assumption seems to hold on very 
large (linear) scales and it is supported by both obser- 



vational evidence (e.g. Baker el al., 1998 and Seaborne 
et al, 1999) and numerical studies (e.g., Blanton et al., 
1999). 

Many authors have used galaxies' peculiar veloci- 
ties and their redshift space positions to estimate the 
value of f3 = Q^'/b, under the hypotheses of linear 
theory and linear biasing. These analyses have been 
typically carried out using two alternative strategies. 
In the so-called density- density comparisons a 3-D ve- 
locity field and a self-consistent mass density field are 
derived from observed radial velocities and compared 
to the galaxy density field measured from large red- 
shift surveys. The typical example is the comparison of 
the mass density field reconstructed with the POTENT 
method (Bertschinger & Dekel 1989, Dekel et al., 1990) 
from the MARK III catalog of galaxy peculiar veloci- 
ties (Willick et al., 1997a) with the galaxy density field 
obtained from the IRAS 1.2 Jy redshift catalog (Sigad 
et al., 1998). The various applications of density- density 
comparisons to a number of datasets have persistently 
led to large estimates of /3, consistent with unity (see 
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Sigad et al, 1998 and references therein). The alter- 
native approach is constituted by the velocity-velocity 
analyses. In this second case the observed galaxy distri- 
bution is used to infer a mass density field from which 
peculiar velocities are obtained and compared to the ob- 
served ones. The velocity-velocity methods have been 
applied to most of the velocity catalogs presently avail- 
able yielding systematically lower values of j3, in the 
range 0.4 — 0.6 (see Zaroubi 2002a, for a summary of 
the most recent results). 

Both density-density and velocity-velocity methods 
have been carefully tested using mock catalogs extracted 
from N-body simulations. They were shown to provide 
an unbiased estimate of the j3 parameter. Yet, when ap- 
plied to the same datasets, the discrepancy in the f3 esti- 
mates turned out to be significantly larger than the ex- 
pected errors. Accounting for mildly nonlinear motions 
(e.g. Sigad et al, 1998 and Willick et al, 1996) or al- 
lowing for possible deviations from a pure linear biasing 
relation consistent with the observational constraints 
(see discussion in Somerville et al., 2001, Branchini 
et a/.2001) does not explain this discrepancy (Berlind, 
Narayanan and Weinberg 2001). Velocity- velocity com- 
parisons are generally regarded as more reliable as they 
require manipulation of the denser and more homo- 
geneous, redshift catalog data. Whereas, the density- 
density comparisons involve manipulation of the noisier 
and sparser velocity data. In any case both classes of 
methods are quite complicated and it is hard to under- 
stand how systematic errors can arise and propagate 
through the analysis. Therefore, it is likely that these 
systematics affect the j3 parameter estimation. 

The purpose of this work is to address, and pos- 
sibly solve, the density- density vs. velocity-velocity di- 
chotomy. We achieve this goal by using the novel Un- 
biased Minimal Variance estimator, recently proposed 
by Zaroubi (2002b). The UMV estimator allows one to 
reconstruct an unbiased cosmological field at any point 
in space from sparse, noisy and incomplete data and to 
map it into a dynamically-related cosmic field (e.g. to go 
from peculiar velocities to overdensities). The UMV is 
applied here to the SEcat catalog of peculiar velocities 
(Zaroubi 2002b) to reconstruct both the mass density 
and peculiar velocity fields. These fields are then com- 
pared with the analogous quantities predicted from the 
distribution of IRAS PSCjz galaxies (Saunders et al, 
2000) of density-density and a velocity-velocity analy- 
ses. In Section H we briefly review the basics of the UMV 
estimator. The SEcat and PSCz catalogs are presented 
in Section ^| Error estimation from mock catalogs is de- 
scribed in § ^. The density and velocity fields obtained 
by applying UMV to SEcat are compared in Section [E| 
with the analogous quantities deduced from the PSCz 
dataset. Finally, in Sections |we discuss the results and 
present our conclusions. 



2 THE UMV METHOD 

The derivation and general properties of the UMV esti- 
mator have been already presented and discussed by 
Zaroubi (2002b). Therefore, we only review its main 



properties. As in the case of Wiener filter (Wiener 1949, 
Zaroubi et al., 1995), the UMV introduces a general 
framework of linear estimation and prediction by mini- 
mizing the variance of the estimator with respect to the 
underlying signal, subject to linear constraints on the 
data. The solution of the minimization problem yields 
the UMV estimator, which was shown to be a very ef- 
fective tool for reconstructing the large scale structure 
of the universe from incomplete, noisy and sparse data 
(Zaroubi 2002b) . 

One of the main drawbacks of the Wiener filter 
is that it suppresses the amplitude of the estimated 
signal. The suppression factor is roughly equal to the 
Signal 2 / (Signal 2 + Noise 2 ) ratio, i.e., in the limit of 
very poor signal-to-noise data, which in the context of 
this work correspond to galaxy peculiar velocities, the 
estimated field approaches zero value. By contrast, the 
UMV estimator has been specifically designed to not al- 
ter the values of the reconstructed field at the locations 
of the data points, thus avoiding spurious suppression 
effects. An unbiased estimate of the reconstructed field 
at any point in space is then obtained by interpolating 
between the data points, according to the correlation 
function assumed a priori. Like the Wiener filter the 
new estimator can be used for dynamical reconstruc- 
tions, i.e., to recover one dynamical field from another 
measured field, e.g., mass over-density from radial pe- 
culiar velocity. These two properties make the UMV es- 
timator a very appealing tool for studying the LSS and 
CMB. 

Here we apply the UMV estimator to reconstruct 
the density and velocity fields in the local universe from 
the radial peculiar velocities of the SEcat galaxies. The 
data points consist of a set of observed radial peculiar 
velocities, u° , measured at positions ri with estimated 
errors ei, assumed to be uncorrelated. The observed ve- 
locities are thus related to the true 3D underlying ve- 
locity field v(r), or to its radial component m , via 

u° = v(r 4 ) ■ Ti + ej = Ui + ti, (2) 

We assume that the peculiar velocity field v(r) and 
the density fluctuation field S(r) are related via lin- 
ear gravitational-instability theory, 8 — /(fi m ) _1 V • v, 
where /(fi m ) ~ f2^ 6 and ftm is the matter mean density 
parameter. Under the assumption of a specific theoreti- 
cal prior for the power spectrum P(k) of the underlying- 
density field, we can write the UMV estimator of the 
3D velocity field as, 

v™ (r) = (v(r)ui)(uittj) _1 w^ (3) 
and the UMV estimator of the density field as, 
8 UMV (r) = ^(pJuiXtn^)- 1 ^. (4) 

Within the framework of linear theory and assum- 
ing that the velocities are drawn from a Gaussian ran- 
dom field, the two-point velocity-velocity and density- 
velocity correlation matrices (bracketed quantities in 
eqs. ^| & H) are readily calculated. Note that the normal- 
ization of the power spectrum drops out of the field es- 
timation. The calculation of these matrices is discussed 
elsewhere (Gorski 1988; Zaroubi et al, 1995,1999). 

The assumption that linear theory is valid on all 
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scales enables us to choose the resolution as well. In par- 
ticular it allows us to use two different smoothing kernels 
for the data and for the recovered fields. In our case no 
smoothing was applied to the radial velocity data while 
we choose to reconstruct the density and velocity fields 
with a finite Gaussian smoothing of radius, R. This 
choice alters the velocity-velocity and density-velocity 
correlation functions that appear in the first bracketed 
terms of the right hand side of eqs. ^ and ^, respectively, 
by introducing a multiplicative term exp[— k 2 R 2 /2] in 
the model power spectrum. 

The amplitude of the reconstructed matter density 
field given in equation ^ is proportional to /(f2 m ) -1 , 
while that of the density field obtained the PSCz galaxy 
distribution is proportional to the biasing parameter, b. 
Therefore, the comparison between these density fields 
will constrain the value of j3 = f(fl m )/b. The velocity 
field reconstruction from the SEcat data set, however, is 
independent of Q m and b. Hence in the velocity-velocity 
comparison {3 enters as a solution of eq. [l] for the PSCz 
velocities, with the matter density given as the ratio 
S(PSCz)/b, where 8(PSCz) is the PSCz density field. 

The error covariance matrix, or variance of the 
residuals, of the reconstructed density and velocity fields 
could be calculated theoretically (Zaroubi 2002b). How- 
ever, in order to give a more complete account of the 
various errors that enter the calculation (e.g., nonlinear 
effects) here we choose to calculate the error from mock 
catalogs (see section [|). 

It is interesting to compare the UMV algorithm 
with other methods of reconstruction used for similar 
purposes. In the context of mass-density reconstruction 
from radial peculiar velocities, comparison with the PO- 
TENT algorithm (Bertschinger & Dekel 1989, Dekel et 
al., 1990) is of special interest. The main assumption 
behind the POTENT algorithm is that the flow field, 
smoothed on large scale, is derived from a velocity po- 
tential. The potential flow assumption is a direct result 
of linear theory but can also be employed in the quasi- 
linear regime until the onset of shell crossing, when an 
extension of eq. [I] applies (Nusser et al, 1991). PO- 
TENT could be viewed as a direct inversion method 
which uses the minimum amount of assumptions, but 
suffers from the problems of direct deconvolution of very 
noisy data. 

In conclusion, the UMV reconstruction can be re- 
garded as a compromise between the POTENT algo- 
rithm, which assumes no regularization but might be 
unstable to the inversion problem of deconvolving highly 
noisy data, and the WF algorithm, which takes into ac- 
count the correlation between the data points and there- 
fore stabilizes the inversion, but constitutes a biased es- 
timator of the underlying field. 



3 THE DATASETS 

The main dataset used here is the SEcat catalog of 
galaxy peculiar velocities which results from the merg- 
ing of 1300 spiral galaxies taken from the SFI catalog 
(Giovanelli et al, 1997a and 1997b and Haynes et al, 
1999a and 1999b) and about 2000 early type galaxies 



from the ENEAR catalog grouped into « 750 indepen- 
dent objects (da Costa et al., 2000). For each object 
the radial velocity and inferred distance, corrected for 
Malmquist bias (Freudling et al, 1995, da Costa et al., 
2000), are provided along with the velocity errors that 
typically amount to ~ 19% of the galaxy distance. 

Merging different peculiar velocity catalogs may re- 
sult in spurious flows and lead to systematic errors in the 
reconstruction procedure. However, several pieces of ev- 
idence indicate that this is not a serious problem for the 
SEcat catalog. First of all, both the SFI and ENEAR 
catalogs are intrinsically homogeneous as they consist 
of uniform data covering most of the sky. Secondly, as 
shown in Bernardi et al. (2002b), the distance relations 
independently calibrated in the two sets are consistent 
with each other and the difference in the Hubble con- 
stant deduced from each catalog is w 5±10 kms _1 /Mpc. 
Finally, different statistical analyses carried out using 
either samples lead to consistent results (e.g., Borgani 
et al, 2001, da Costa et al, 2001, Nusser et al., 2001, 
Zaroubi et al., 2001 and references therein). Some of the 
main characteristics of these samples are summarized in 
Figure |l| The sky distribution (upper panels) , redshifts 
(central panels) and peculiar velocities (lower panels) of 
the SFI and ENEAR galaxies are quite similar, espe- 
cially when accounting for the expected tighter spatial 
correlation and higher peculiar velocities of the ENEAR 
early type galaxies. Further evidence for the consistency 
of the two catalogs will be provided in section ^ (e.g., 
see Figures ^ and . 

The substantial number of galaxies and the large 
sky coverage (the unobserved region is given by a zone 
of avoidance of about 15° about the Galactic-plane) al- 
low a dense and uniform sampling of the velocity field. 
Moreover, since SEcat contains both early and late type 
galaxies, we can sample both high and low density re- 
gions and therefore minimize possible biases that might 
have affected other analyses based on a single popula- 
tion of objects. 

To estimate the /3 parameter one needs to com- 
pare the densities and/or velocities reconstructed from 
a radial velocity survey with those recovered from a di- 
rect probe of the density field, namely a redshift sur- 
vey. Here we use the models obtained by Branchini et 
a/.(1999) from the distribution of IRAS PSCz galaxies 
under the assumptions of linear biasing and linear the- 
ory. The PSCz redshift survey (Saunders et al, 2000) 
provides the angular positions and redshifts of ~ 15000 
IRAS galaxies distributed over almost all sky (the zone 
of avoidance is about 8° ) with a median redshift of 8500 
kms _1 and is therefore suitable for modeling the density 
and velocity fields within the same region in which the 
UMV reconstruction is performed. The typical error as- 
sociated with the PSCz density and velocity models are 
significantly smaller than those reconstructed from the 
SEcat velocities and therefore will be ignored in the sub- 
sequent analysis. 
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Figure 1. Upper panels: The projected distribution of the SFI (left panel) and ENEAR (right panel) galaxies in galactic 
coordinates. Crosses indicate positive peculiar velocities, open circles negative; and the size of the symbols is proportional to 
the amplitude of their peculiar velocity. Middle panels: The redshift distribution of the SFI (left panel) and the ENEAR (right 
panel) galaxies. Lower panels: The peculiar velocity distribution of the the SFI (left panel) and the ENEAR (right panel) 
galaxies. 



© 0000 RAS, MNRAS 000, 000-000 



SEcat vs. PSCz 



Underlying density 




Mock without errors 



oe 

SGX (Mpc/h) 



Mock with errors 



s 





-50 <E 50 

SGX (Mpc/h) 



Mock with errors 




-50 OB 

SGX (Mpc/h) 



50 



-60 



oe 

SCX (Mpc/h) 



50 




1 

<5(0riginal) 



Figure 2. Comparison of the original and reconstructed G12 smoothed density maps on the mock Supcrgalactic plane. In all 
panels the solid and dashed line contours denote positive and negative densities respectively. The bold-solid line denotes the 
zero level density. Contour spacing is 0.1. The very thick solid contour marks the area within which the reconstruction errors 
are less than 0.2. The upper left hand panel shows the underlying density field of the simulation. The degradation of the original 
density map towards the edges is spurious and due to the finite size of the smoothing radius. The upper right hand panel shows 
the reconstructed density from the SEcat mock catalog before adding errors to the distances and velocities. The other three 
maps show reconstruction from SEcat mock catalogs with realistic noise. The lower right-hand panel shows a typical scatter 
plot of the density within the comparison region of the original vs. the reconstructed density from one of the mock catalogs 
(with noise); the points used in the scatter plot are 1/10 randomly sampled from the grid points with estimated reconstruction 
errors less than 0.2. 



4 MOCK CATALOGS AND ERROR 
ESTIMATES 

To test the performance of the method when applied 
to the SEcat catalog, we construct mock catalogs based 



on the "Constrained Realization GIF" simulation car- 
ried out by Mathis et al. (2001). This simulation starts 
from initial conditions with a smoothed linear density 
field which matches that derived from the IRAS 1.2 Jy 
galaxy survey and tracks the formation and evolution 
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Figure 3. The distribution of the value of f) relative to the real one as estimated from 100 mock SEcat catalogs in the density- 
density comparison (left panel) and velocity- velocity comparison (right panel). The value of the measured /3 for the catalog 
without noise is (3 = 0.99 and 1.01 for the density and velocity reconstruction respectively. 



of all dark matter halos more massive than 10 11 solar 
masses out to a distance of 8000 km s^ 1 from the Milky 
Way up to the current epoch. Galaxies in the original 
mock catalog are sampled from the N-body simulation 
at the same position of the real galaxies. Realistic mock 
catalogs are obtained by assigning errors to the position 
and velocity of mock galaxies consistent with observa- 
tions. We construct 100 mock catalogs that differ only 
in the realization of errors added to the position and 
velocity of the galaxies. Measured peculiar velocities in 
the mock catalogs are then obtained after performing a 
Malmquist bias correction according to the recipe given 
by Willick et al. (1997a). 

The density and 3D velocity fields within a region 
of 60 h~ x Mpc , smoothed with a Gaussian filter of ra- 
dius of 12 h' 1 Mpc (G12 hereafter), were reconstructed 
from each of the 100 mock catalogs. Both N-body and 
reconstructed density and velocity fields were specified 
on a regular grid with a mesh size of 2.5 Mpc and re- 
constructed radial velocities at the actual location of the 
data points were interpolated from the grid. Errors in 
the reconstructed densities and peculiar radial velocities 
were estimated from the 100 Monte Carlo realizations. 

Note that we do not use the peculiar velocities of 
the mock galaxy catalogs of Mathis et al. (2001). In- 
stead, the peculiar velocities of our mock galaxies are 
obtained from those of of the dark matter particles. This 
implies neglecting the so called velocity bias which is ex- 
pected to be small on the large smoothing scale involved 
in our reconstruction (see e.g., Carlberg 1994). 

Figure shows the quality of the reconstruction 
from mock SEcat catalogs. The comparison of the origi- 
nal underlying density with the one reconstructed from 
ideal mock catalogs (i.e. with no velocity errors) in a 



region in which the estimated reconstruction errors are 
< 0.2 shows that their relation is well described by a 
linear function with a slope of 0.99 (see below). This 
demonstrates that the sampling density of the data is 
sufficiently high. Note that the density of the original 
catalogs close to the boundaries is suppressed. This ef- 
fect is due to finite smoothing length and does not af- 
fect our analyses which are carried out to distances well 
within the edges of the mock catalogs. The other three 
maps shown in Figure 2 were randomly chosen from the 
100 mock catalogs reconstructions. The similarity of the 
reconstructions are quite evident, especially within the 
region of small reconstruction errors. The scatter plot 
quantifies this similarity for one of the mock catalogs. 
An inspection of the radial peculiar velocity reconstruc- 
tion shows a similar quality of results. 

Left panel of Figure ^ shows /3 as estimated by ap- 
plying the UMV reconstruction algorithm to the mock 
catalogs and comparing the results with the true density 
and velocity fields using the following x 2 statistic, 

2 _ l_ ST [^(Mock) - /3(5i (Original) + AS)] 2 
X N 2^ a \ ' 1 ' 

cr s <0.2 

where Si (Mock) is the density as reconstructed from the 
mock SEcat catalog; <5; (Original) is the original density, 
and as is the rms difference between the reconstructed 
densities from the mock catalogs and the original den- 
sity at the same point in space. The free parameters 
here are /? and the offset between the two fields, AS. 
The latter is introduced to account for the uncertainty 
in denning the mean density of the sample. This off- 
set is expected to be zero in the mock catalog analyses 
and can only be non trivial when comparing the den- 
sity fields obtained from two different catalogs (such as 
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the true SEcat and PSCz catalogs). The density-density 
comparison is carried out over N points, randomly se- 
lected from the 10 x N at which the estimated errors 
are less than 0.2. The expected values of /3 and AS are 
unity and zero respectively. The left panel of Figure |^ 
shows that the estimated (3 is unbiased and has a mean 
and variance of 1.01 and 0.1, respectively. The estimated 
value of A<5 is —0.03 and has a variance of 0.04. 

Naturally, most of the grid points used in the com- 
parison - even after diluting their number by a factor of 
10 - do not have statistically independent errors. There- 
fore, the statistic used in eq. [| is not optimal. In other 
words, had the data points used in eq. |E] been indepen- 
dent then it would indeed be sufficient to use the likeli- 
hood contours obtained after minimizing the x 2 statistic 
to determine the uncertainty of the results. On the other 
hand, had these data point been totally dependent with 
about only one degree-of-freedom then the error ob- 
tained from the likelihood analysis would be a gross un- 
derestimate unless the cosmic-variance like errors over 
the whole sample were also estimated and taken into 
account. Normally, the best way to go about estimat- 
ing the scatter in the results that takes into account 
the partial dependency of the data points and their lim- 
ited number is to perform a very large number of Monte 
Carlo simulations - of the order of 10 2 -10 3 - to ensure 
an accuracy of a few percent. This however is very time 
consuming and not feasible with available computer re- 
sources. Therefore, we choose to assume that the errors 
are totally independent and use eq. [| in order to esti- 
mate the most likely /3 and AS values. The uncertainties 
are then determined by adding in quadrature the like- 
lihood errors and the uncertainties obtained from the 
scatter in the value of /3, obtained from the 100 mock 
catalogs. In our sample, these two sources of errors are 
not totally independent. However, by treating them as 
such allows one to obtain a conservative upper limit on 
the errors. 

In order to validate our approach it is important to 
estimate the number of degrees-of-freedom in the sam- 
ple, Mdof- In the case of no noise correlations except 
those introduced by the G12 smoothing, Mdof ~ 23 in 
a sphere of radius bOh^ 1 Mpc , which represents the ra- 
tio of the total volume to the effective volume of the 
G12 filter. In the case at hand, however, the calcula- 
tion is more subtle and involves numerical estimation of 
the noise correlation matrix from the 100 Monte-Carlo 
simulations at each grid point with error less than 0.2 
and finding its eigenvalues. Then Mdof is identified with 
the number of significant eigenvalues of this matrix (see 
Zaroubi et at, 1995 for the treatment of a similar prob- 
lem). Specifically, Mdof is the number of the highest 
eigenvalues that account for 95% of the variance, found 
here to be about 20 (see Figure ^]). The eigenvalues of 
most of the remaining eigenmodes drop by orders of 
magnitude. This number, 20, reflects the number of ef- 



The density reconstruction is performed hereafter normal- 
izing with the correct value of f(Q m = 0.3) therefore the 
comparison is expected to yield a bias parameter of unity. 
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Figure 4. The sorted spectrum of the eigenvalues of the 
noise correlation matrix. The modes that account for 95% of 
the variance are the first 20. 



fective smoothing volumes and the additional correla- 
tion introduced by the UMV filter. 

The same strategy is adopted for the velocity- 
velocity comparison. In this case the statistic of choice 
is, 



N d , 



E 

DataPoints 



(ui( M °ck) - 0Ui (Original) - AH n 



where u(Mock) and w(Original) are the radial veloci- 
ties of the mock and original data respectively. Ndata is 
the number of data points, AH a is the offset in Hub- 
ble constant, is the distance of the data point i and 
<r v is the uncertainty in the radial velocity as obtained 
from the mock catalogs. The velocity-velocity compari- 
son is carried out over the radial velocities at the loca- 
tion of the data points. The argument regarding the \ 2 
statistic used here is identical to the one discussed for 
the density-density comparison. Right panel of Figure |^ 
shows that the mean and variance of f5 in this case are 
0.97 and 0.12, respectively. The offset in Hubble con- 
stant is 0.2 ± 0.5kms _1 /Mpc. The value of /3 estimated 
from the noise- free mock catalog is 1.03. From the noise 
correlation matrix we estimate Mdof ~ 17. The lower 
value obtained here reflects the longer range of velocity 
correlation. 

Willick and Strauss (1997b) have suggested in their 
VELMOD method to estimate /3 from velocity data 
by minimizing residuals and their correlations simul- 
tanously with minimizing the likelihood, therefore, al- 
lowing the calibration of the distance estimator while 
evaluating (3. Whereas in our analysis we assume that 
the distance indicator has been calibrated indepen- 
dently. 

The analysis of the mock catalogs shows that the 
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Figure 5. The G12-smoothcd, overdensity field on the Supergalactic plane reconstructed from SEcat, ENEAR and SFI catalogs. 
All reconstructions assume a Standard CDM power spectrum apart from the map shown in the upper right panel for which we 
have assumed an Q m = 0.3 ACDM power spectrum. The very thick solid contour marks the area within which the reconstruction 
errors from the SEcat catalog are less than 0.2. In all panels the solid and dashed line contours denote positive and negative 
densities respectively. The bold-solid line denote the zero level density. Contour spacing is 0.1. 



UMV returns an unbiased estimate of the underlying 
density and velocity fields from which 13 can be deter- 
mined with an accuracy of 10%. This error is purely 
random and takes into account the scatter of the best 
fit /3 value but does not include the scatter in the esti- 
mation of each individual value of f3 which will be added 
later. 

This analysis could in principle be extended to 
smaller smoothing kernels. Besides the validity of lin- 
ear theory required in this analysis, the density of the 
sky coverage might be also an important factor since 
an insufficient coverage might result in the dominance 



if noise in most of the sampling volume. An extension 
of the analysis to a 900 km s -1 smoothing kernel (here- 
after, G9) shows that the results of the analysis are con- 
sistent with the G12 comparison but the uncertainties 
are much higher especially in the density-density com- 
parison. Therefore, the smoothing kernel used in the rest 
of the paper is G12. 
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Figure 6. A quantitative comparison between the SFI and 
ENEAR G12-smoothed. reconstructed densities. The densi- 
ties of both SFI and ENEAR were reconstructed on a grid 
with mesh size of 2.5 h~ 1 Mpc The densities shown in the 
scatter plot are from gridpoints randomly selected with a 
rate 1/10 from those with reconstruction errors less than 
0.2. The size of the symbols is inversely proportional to their 
errors. The solid line with a slope of unity has been drawn to 
guide the eye. The agreement between the two reconstruc- 
tions is very good except for a small number of points with 
<5(SFI) ss -0.5 and <5(ENEAR) ss 0. 

5 THE FIELD-FIELD COMPARISON 

Here we apply the UMV estimator to the true SEcat 
catalog to obtain the G12-smoothed density and veloc- 
ity fields assuming a power spectrum of a flat CDM 
model with h = 0.5, Q m = 1.0 as a prior. These cos- 
mological parameters determine the shape of the power 
spectrum only and are not used anywhere else in the 
analysis. Unlike the outcome of the Wiener Filter es- 
timator, the UMV-reconstructed density and velocity 
fields, for dense sampling of the sky, are unbiased and 
therefore can be used for quantitative comparisons. 

Before proceeding further in analysing the SEcat 
catalog, it is useful to compare the reconstruction we 
obtain from the SEcat peculiar velocities with those ob- 
tained from the SFI and ENEAR separately in order to 
support our claim regarding the consistency the latter 
two. Figures ^, |(| and ^ show the similarity between the 
density and velocity reconstruction of the three cata- 
logs. We also note that, as demonstrated in the same 
figures, changing the power spectrum prior has a very 
small effect on the results. 

The reconstruction of the PSCz model density and 
velocity fields, performed according to Branchini et 
a/.(1999), requires an input value for j3; here we use 
(3 = 0.5. However, the dependence of the reconstructed 
density field on the input (3 is very weak while the 
model peculiar velocity roughly scales linearly with (3. 



This dependence results in a systematic error of about 
2% on the final result. Model PSCz velocities were ob- 
tained at the reconstructed positions of PSCz galax- 
ies in real space. Both masses and velocities were then 
smoothed with a G12 filter to obtain the PSCz density 
and velocity fields on a regular grid with a mesh size 
of 2.5 h' 1 Mpc Finally, G12-smoothed PSCz velocities 
have been interpolated at the positions of SEcat galax- 
ies. Here we compare the SEcat fields with those of the 
PSCz, all smoothed with a G12 filter. The errors on the 
model density and velocity fields are much smaller than 
the SEcat ones (see Branchini et al, 1999) and will be 
ignored in the following comparisons. 

5.1 The density-density comparison 

The left panels of Figure ^ show the G12 smoothed 
UMV reconstructed density field map in three planes 
at different Supergalactic Z (the central plane refers to 
Z=0, i.e. the Supergalactic plane) obtained from the 
SEcat peculiar velocities within a box of 160 h~ x Mpc 
aside, centered around the Local Group position. The 
main features of our local universe are easily identified 
in the UMV map on the Supergalactic plane , including 
the Great Attractor on the left and the Perseus-Pisces 
supercluster in the lower right. There is also a hint of 
the Coma cluster, which lies just outside the sample, in 
the upper part on the map. Similar features also char- 
acterize the PSCz density map shown on the right-hand 
panel of Figure ^| This map, obtained by Branchini et 
a/.(1999), has the same smoothing (G12) and shows the 
same region of the universe. 

A quantitative density-density comparison is car- 
ried out using the following \ 2 statistic: 

2_}_ ST (SEcat) - /3(^(PSCz) + AS)] 2 
X N 2^ al ' K ) 

a s <0.2 

where Si(SEcat) and <5i(PSCz) are the SEcat and PSCz 
densities respectively, AS is the offset in the mean den- 
sity and as are the density reconstruction error at each 
point as estimated from Monte Carlo realizations of 
mock-SEcat catalogs. The best fit f3 and A<5 parame- 
ters are those that minimize Eq. |^. Here again the sum 
is over gridpoints randomly sampled at a rate of 1/10 
from those that have errors less than 0.2. As in the mock 
catalogs analysis, here Ndof ~ 20. 

In previous density-density comparisons {e.g., 
Sigad et al., 1998) the authors chose to minimize the 
X 2 with respect to 1/(3 instead of (3 directly. Since the 
main source of errors in our analysis are the uncertain- 
ties in the measured galaxy velocities, adopting a direct 
or inverse f3 minimization of the x 2 statistic does not 
make much difference. Here we choose to minimize with 
respect to (3. 

Each point in the scatter plot displayed in the lower 
panel of Figure ^| shows the comparison between the 
SEcat and PSCz overdensities, measured at the same 
locations. The comparison is restricted to the 1/10 ran- 
domly chosen points with as < 0.2. The slope of the 
solid line gives (3 = 0.57. 

A zero-point offset AS = 0.18 is also detected. We 
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Figure 7. The top left panel compares the G12-smoothed peculiar radial velocities reconstructed from the SFI catalog versus 
those reconstructed from the SEcat catalog, both measured at the location of the SFI data. The solid line represents the best 
linear fit to the scatter plot. The top-right panel compares ENEAR and SEcat reconstructed velocities at the positions of 
the ENEAR galaxies. The two bottom panels refer to the SEcat catalog only and compare velocities reconstructed assuming 
a ACDM model (Y axis) with those reconstructed from the standard CDM model (X axis). Peculiar velocities shown in the 
bottom-left panel refer to points with estimated error less than 100/fcms. All points are considered in the bottom-right panel. 
In all plots the size of the symbols is inversely proportional to their errors. 



interpret it as a mismatch in the average density in the 
two samples. The mismatch in the mean fields is caused 
by the PSCz density field which was found to be system- 
atically larger than the IRAS 1.2Jy density field within 
a 60/i _1 Mpc sphere (Teodoro et al, 2000). This mis- 
match arises due to the incompleteness of the PSCz 
catalog at low fluxes. 

The upper panel in Figure ^| shows the 1, 2 and 3- 
a likelihood contours in the /3-A5 plane obtained from 



Eq. 0. The marginalization of this distribution with re- 
spect to AS gives the error estimate on the value of /3 
which is of the order of 0.08 (~ 15%). Adding to this 
error the error estimated from the distribution of bias 
shown in Figure ta one obtains f3 = 0.57+° ^3. 

To estimate the goodness-of-fit of the parameters 
obtained from the \ 2 analysis, we calculate the distri- 
bution of the residuals, 
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Figure 8. The panels on the left show the maps of the G12-smoothed, density fluctuations UMV-reconstructed from the 
SEcat catalog. The central panels show the same field reconstructed from the distribution of PSCz galaxies. The maps of 
density residuals, computed for /3 = 0.57 A<5 = &SEcat ~ SpsCz are shown on the right-hand panels. The central maps 
show the density fields on the Supergalactic planes. The maps on the upper and lower panels show the density fields at 
Supergalactic Z= ±20 _1 Mpc, respectively. The very thick contour marks the boundaries of the volume within which the 
estimated reconstructed error is < 0.2; the volume used for comparisons. 



<5(SEcat) - Pml(S(PSCz) - AS) 



where /3ml is the best fit j3 parameter and as is the 
error on the UMV estimated density field. If the model 
correctly describes the data, this distribution should be 
Gaussian with a rms of unity. The histogram of £ is 



shown in Figure [10| along with the best fitting Gaussian 
distribution (dashed line), whose rms is almost unity. 

To check the robustness of the result we have re- 
peated the density-density comparison by drawing three 
additional samples defined at different error thresholds 
and re-estimating the values of ft. Adopting the thresh- 
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Figure 10. The histogram represents the distribution of 
the normalized residuals (£). The dashed line shows its best 
Gaussian fit with a rms value of <r = 1.03, as indicated in the 
plot. 
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-1 



(5(PSCz) 

Figure 9. The upper panel shows the 1, 2 and 3 a likelihood 
contours from the 5 — 5 comparison in the /3 — A<5 plane. The 
scatter plot in the bottom panel compares the G12-smoothed 
UMV-reconstructed SEcat density field to the PSCz density 
fluctuations, both measured at grid points 1/10 randomly 
selected from those with reconstruction errors ug < 0.2 The 
size of the symbols is inversely proportional to their errors. 



olds as = 0.3, 0.4 and 0.5 we obtained /3 
and 0.52, respectively. 



0.54,0.53 



5.2 The velocity-velocity comparison 

The UMV- velocity reconstruction procedure is different 
from the one used for the density. SEcat radial veloci- 
ties are calculated at the location of the data points from 
the reconstructed 3D velcoty which has been homoge- 
neously smoothed with a G12 filter through the UMV 
operation. Here we compare these velocities to the G12- 
smoothed model PSCz velocity sampled at the same lo- 
cations. Therefore, unlike for the density-density case, 
the number of points used in the velocity-velocity com- 
parison is determined directly by the number of points 
in the SEcat catalog. This comparison is carried out 
using the following % 2 statistic 

\2 



1 

N 



(m (SEcat) - f3u t {PSCz) - AH a n 



DataPoints 



where u denotes the G12-smoothed radial velocities, 
AH a is a local perturbation to the Hubble constant, 
7-j is the radial distance of the point i and a v is the er- 
ror in the UMV reconstruction. The resulting velocity- 
velocity scatter-plot is shown in the lower panel of Fig- 
ure [ll] along with the best fitting line. Here the slope of 
the line constitutes an estimate of j3 = 0.51. The zero- 
point mismatch, AH a represents a spurious "breathing- 
mode" which is to be expected given the average den- 
sity mismatch found in the density-density comparison. 
A perturbation of AH a = 1.5km s -1 /Mpc was found 
and its associated spurious radial motion, AH a r, was 
subtracted from the PSCzneculiar velocities shown in 
the lower panel of Figure [ill. The two zero-points are 
found here to be consistent with the prediction of the 
linear theory relation: 
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Figure 11. Upper Panel: The 1,2, and 3 a likelihood con- 
tours in the — AH o plane from the SEcat vs. PSCz velocity- 
velocity comparison. In the bottom panel the G12-smoothed, 
reconstructed SEcat and PSCz are compared at the locations 
of the SEcat data points. The size of the symbols is inversely 
proportional to the reconstruction errors 



ft ' 6 

AH = -J^AS(< r)H , 



(10) 



where A<5(< r) is the mean-density mismatch within a 
radius r. 

The upper panel of Figure [ll] shows the likelihood 
contours in the f3 — AH plane obtained by calculat- 
ing the x distribution given in Eq. ^. From the iso- 
probability contours of 1,2 and 3-<r levels shown in the 
upper panel of Figure [n] we obtain that /3 = 0.51 ±0.06, 



o 

u 




-6 -4 -2 2 4 6 
[v(SEcat)-/?v(PSCz) -AH ]/ct v 

Figure 12. The distribution of the residuals from the veloc- 
ity velocity comparison. The dashed line shows a Gaussian 
fit with cr = 1.06. 



fully consistent with the estimate of /3 from the density- 
density comparison. The uncertainty here is a combina- 
tion of the errors estimated from Figures ]§] and [ll] 

The analysis of the velocity-velocity residuals is 
performed similarly to those of the density-density resid- 
uals. The resulting distribution of the residuals is shown 
in Figure ^ as a histogram. Here again the rms value 
found for the best fit Gaussian distribution (dashed line) 
is very close to unity indicating the adequacy of the 
PSCz velocity model with f3 = 0.51. The slight excess 
in the positive tail of the histogram is due to presence of 
few outliers,easily identified below the best fitting line 
in the lower panel of Figure [ll]. The size of those points 
clearly indicates that they have very large errors, which 
would exclude them from any reasonable noise data- 
cut analysis. Moreover, their number is very small and 
therefore their influence on the final result is negligible. 



6 SUMMARY AND DISCUSSION 

In this paper we have applied the new UMV estimator 
to recover the density and velocity fields in the local uni- 
verse from the SEcat catalog of galaxy peculiar veloci- 
ties. In order to obtain the so called (3 parameter, these 
fields have been compared with those modeled from the 
spatial distribution of IRAS PSCz galaxies assuming 
linear theory and biasing. Previous estimates of f3 from 
density-density comparisons, mainly based on the PO- 
TENT algorithm (Bertschinger & Dekel 1989, Dekel et 
al, 1990), have yielded a large value (/3 » 1 c/., Sigad et 
al, 1998), inconsistent with the smaller values (« 0.5) 
independently obtained from all recent velocity-velocity 
VELMOD (Willick et al, 1996, Willick & Strauss 1998, 
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and Branchini et al, 2001) and ITF (da Costa et al, 
1998 and Nusser et al, 2000) comparisons. 

For the first time the UMV method provides a 
common methodological framework in which to perform 
velocity-velocity and density-density comparisons. The 
velocity-velocity comparison yields a value of (3 con- 
sistent with that measured in the VELMOD and ITF 
analyses. However, the value of the same parameter ob- 
tained from our density-density comparison is signifi- 
cantly smaller than those obtained from the POTENT 
analyses (c/., Sigad et al, 1998). The f3 parameters from 
both v — v and 8 — 8 comparisons presented here are in 
agreement, yielding a /3 w 0.55 with an estimated error 
of the order of 0.1. 

In contrast with the POTENT algorithm, the new 
UMV method reconstructs the density field from pecu- 
liar velocities while taking into account their underlying 
correlation properties. The regularization aspect of the 
UMV estimator significantly improves the stability of 
the inversion, which is especially important given the 
low signal-to-noise ratio of peculiar velocity data. The 
regularization obtained by this method is very similar to 
the one provided by the Wiener filter method (Zaroubi 
et al, 1999). However, the UMV has been designed to 
provide an unbiased estimator of the underlying signal, 
a property that the Wiener filtering method lacks. These 
two aspects make the UMV estimator a very useful tool 
for reconstruction from peculiar velocity data. 

In our error analysis we have shown that for the 
best fit value of (5 the density and velocity residuals 
are normally distributed. This indicates that the PSCz 
density and velocity fields constitute an adequate model 
for those reconstructed with the UMV estimator. The 
fields only differ by a monopole term, corresponding to 
a mismatch in the the mean density within 60 /i -1 Mpc 
which is caused by the known incompleteness of the 
PSCz catalog at faint fluxes (Teodoro et al, 2000). This 
also implies that the effect of the nonlinear dynamics 
and amount of nonlinear and stochastic biasing on the 
scales involved in our analysis is negligible relative to 
the measured peculiar velocity errors. 

The results presented in this paper are quite en- 
couraging since for the first time the two ways of esti- 
mating the value of the (3 parameter give a consistent 
result. Our results also suggest that the UMV estimator 
is a promising tool for the problem of reconstructing the 
dynamical fields from the observed radial peculiar veloc- 
ities and therefore could be applied to other datasets. In 
particular, to reconstruct the large scale structure from 
the incoming large and uniform surveys that will provide 
both the spatial distribution and peculiar velocities of 
extra-galactic objects e.g., SDSS and large cluster sur- 
veys with kinematic Sunyaev-Zel'dovich measurements. 

Our present density- density comparison results are 
in marked contrast to those obtained by earlier work, 
including those from the recent analysis of the Mark III 
catalog using the POTENT method {e.g., Sigad et al, 
1998), and raises the question on the origin of this dis- 
crepancy. The Mark III catalog, as shown for example 
by Davis et al. (1996) and more recently by Courteau et 
al. (2000), suffers from systematic calibration errors that 
would cause a systematic error in the estimation of (3. 



However, these errors are not expected to overestimate 
the value of f3 by more than a factor of two. An applica- 
tion of the UMV method to the Mark III catalog shows 
that the obtained values of (3 are somewhat higher than 
those obtained from the SEcat catalogs by 0.1 — 0.2. 
Moreover, the v — i>-like VELMOD analysis yield con- 
sistent values of (3 when applied to Mark III and SFI 
datasets (Willick et al., 1997b, Willick and Strauss 1998, 
Branchini et al., 2001). Based on these arguments, one 
could speculate that the most likely explanation to the 
inconsistent results is a conspiracy of both the system- 
atics errors in Mark III and some noise-driven inversion 
instability in the POTENT reconstructions. 
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